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■ We investigate Isaacson's high-frequency gravitational waves which propagate in some 

relevant cosmological models, in particular the FRW spacetimes. Their time evolution in 
^ |. Fourier space is explicitly obtained for various metric forms of (anti-)de Sitter universe. 

D ' Behaviour of high-frequency waves in the anisotropic Kasner spacetime is also described. 
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^ ■ 1 Introduction 

. During the last years a growing attention has been paid to analysis of gravita- 

tional-wave perturbations propagating on cosmological backgrounds. The motiva- 
tion stems not only from theoretical considerations concerning the influence of 
' primordial gravitational waves on the cosmological evolution, in particular on the 

. CMB anisotropics, but also from the emerging field of gravitational- wave astron- 

' omy which stimulates studies on stochastic backgrounds of gravitational waves (see, 

Q^i e.g., [1-10] for number of references). 

^ I Many years ago Isaacson introduced and developed an interesting gauge invari- 

bX). ant perturbation approach that enables one to study properties of high-frequency 

gravitational waves, together with their influence on the background in which they 
propagate [11,12]. The complete Isaacson's formalism thus consistently combines 
' the equation for high-frequency waves with the equation describing back-reaction 

of the waves on the background. Recently we investigated this effect using the 
WKB form of the gravitational-wave perturbations [13] and we presented relations 
between these solutions and some exact radiative cosmological models. 

It is rather surprising that the Isaacson wave equation itself did not attract 
much attention. To our knowledge it was solved (in fact, by several authors [11-16]) 
basically in case of the Robinson- Trautman spacetimes, the Vaidya solution in par- 
ticular. However, such spacetimes are usually considered as a model of spherical 
gravitational waves generated by a "local source" rather than a cosmological model. 
The purpose of our present work is to solve the Isaacson wave equation describ- 
ing propagation of high-frequency waves on "genuine" cosmological backgrounds 
(without, however, considering the back- reaction). 

After summarizing the Isaacson formalism (section 2) we apply it to some 
fundamental cosmological models (section 3). General family of the Friedmann- 
Robertson- Walker spacetimes is studied in subsection 3.1. By decomposing the 
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high-frequency perturbations into tensor harmonics we derive the evolution equa- 
tion for various permitted gravitational waves. In particular, illustrative explicit 
solutions are presented for the de Sitter and anti-de Sitter universes using var- 
ious metric forms of these spacetimes. In subsection 3.2 we investigate possible 
high-frequency gravitational waves in the Kasner universe which is a classic repre- 
sentative of anisotropic but homogeneous cosmological models. 



2 The Isaacson formalism 



Isaacson's perturbation method [11] is based on decomposition of the spacetime 
metric g^u into the background metric 7^^!/ and its perturbation hfj,^, 

where, in a suitable coordinate system, 7^^ — 0(1) and h^i, — 0(e). By defini- 
tion, / — 0(e") if there exists a constant C > such that |/| < Ce" as e ^ 0. 
The assumption /i^^ = 0{e) thus does not automatically imply that /i^^ ~ e. The 
spectrum of possible waves is not a priori restricted, it is only required that their 
amplitudes fall to zero at least hnearly with e, i.e. |ft.^i/(e)| < Ce. 

Physically, the non-negative dimensionless parameter e is the ratio of a typical 
wavelength A of gravitational waves and the scale L on which the background 
curvature changes significantly. Isaacson's high-frequency approximation thus arises 
when A ^ L, i.e. e ^ 1. Consequently, 

7m. = 0(1) , h^, - 0(e) , 

7M.,a = 0(1) , V,o = 0(1) , (2) 

liiu,a0 = 0(1) , h^^^al3 = 0(e"^) , 

and one may expand the Ricci tensor in powers of perturbations h as 

R,. (g) = <J + <j + Ri'} + ■■■ , (3) 

where 

/?(!) ^ L^P'^ (h j^h —h \ (A\ 

^'fiv — 2 ' \"'T/j.;i>p ~ '^TU-.iAp '^pT-.fiiy '^p^u^pr} 1 y*) 

IhP-^ h I /jP^f/j 4-/7 — /7 _ /7 ^ 

2'^ -^f'^pr^fj, \ \'^Tp-.fj,iy I "'fiiy;Tp '^Tp\yp "'Tiy.p,pJ 



i?(2) ^ 1 



The semicolons denote covariant differentiation with respect to the background 
metric 7^^ , which is also used to raise or lower all indices. Considering the relations 
(|2Jl, the orders of the terms (0J) are 

4°^=0(1), <j=0(e-i), i?(l)=0(l). (5) 
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In the Isaacson high-frequency approximation e <C 1 the dominant term is R^}} 
yielding the wave equation 

{hrn'^up ^ru-^^p hp'r-^pu h^j^-^pr^ — (6) 

for the perturbations hpi, on the curved background (considering the vacuum 
complete metric gp^, or the case of a stress-energy tensor which does not contain 
derivatives of the metric, as shown in [13]). The two terms of the next order 0(1), 
namely R!^p} and r!]}} , can be used to provide the equation for the background (non- 
vacuum) metric [11,13] which involves the essential influence of the high-frequency 
gravitational waves on the background. Of course, to obtain a consistent solution, 
one has to solve both the wave equation and the equation for the background 
simultaneously. We analyzed this problem explicitly in [13] for backgrounds with 
privileged non-twisting null directions, after simplifying the equations by the WKB 
approximation. 

In this article, however, we wish to concentrate on the principal wave equation 
(|HJl. Our aim is to obtain evolution of high-frequency gravitational radiation prop- 
agating in an arbitrary direction in various cosmological models. If we impose the 
TT gauge conditions 

V'' = 0: (7) 
h^t. = 0, (8) 

the equation ((HJ reduces to the following wave equation 

- 24°j^^ h^" - h% - <j = . (9) 

Consistency of the gauge conditions with the wave equation ^ was demonstrated 
in [11]. The conditions © still do not completely exhaust the gauge freedom, 
and we can thus prescribe the following additional condition 

hpo = (10) 

to further simplify calculations. 

3 High-frequency waves in cosmological models 

Our aim is to investigate solutions of the wave equation lO, subject to the 
gauge conditions |2l and (O, for some relevant cosmological models, namely the 
Friedmann-Robertson- Walker (FRW), (anti-)de Sitter, and anisotropic Kasner 
universes. 

3.1 Waves in the FRW models 

A general FRW metric can be written in the standard form 

ds^ = a^{r]){-drf + dx'dx^) , (11) 
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where rj is the conformal time and ■^7^ is the metric of homogeneous and isotropic 
three-space with the uniform spatial curvature K. Usual choice of the spatial coor- 
dinates is 

dx'dx^ = dx^ + .g^(x)(d6'2 + sin^ 9d<j)^) , (12) 

where g — sinx for K = 1, g — sinhx for K = —1 and g = x for K = 0. Let 7^/3 
and ; denote the FRW metric (|ll(l and the corresponding covariant derivative, 
respectively. The covariant derivative with respect to the metric '^7^ will be denoted 
by |. 

We consider a specific form of the metric perturbations (which has been already 
used in literature, see e.g. [1]): the time and spatial dependence of the perturbations 
are separated in the following way, 

= fiv) Q^Lu , (13) 

where satisfies Q^o = 0, in accordance with the gauge condition (|l()|l . The 
spatial components of Qij form a traceless divergenceless tensor (thus guaranteeing 
that hfj^i, satisfies the gauge conditions iQ, ©) which is a solution of 

+ =0 . (14) 

Such functions Qij are called tensor harmonics, and equation (|14(l is a generalized 
Helmholtz equation with k representing the wave number which sets the scale of 
the perturbations relative to the background coordinates. The expansion of pertur- 
bations into the tensor harmonics was investigated from the mathematical point of 
view, e.g., in [17]. 

To simplify the form of the curvature terms in the wave equation one can 
use the well-known decomposition of the Riemann tensor [18] 

RafivIS = Ccrfj.up + la[uRl3]ti - lti[uRp]a ~ -^la[ylp\nR , (15) 

where Ca^vp is the Weyl tensor. For the FRW metric pi(l the Weyl tensor vanishes, 
and one obtains 

i?y = a^'^{aa + + 2X0^)7^ , R ^ 60"'"^ (a + Ka) , 

where R = ^^'^R^i^- From H15|l . (jSl, H10() . and the fact that the spatial part Rij of 
the Ricci tensor is a multiple of 7ij , we thus derive 

Ra^.phi''' = \ (r% + R\ ~ \r) V (16) 

(no summation over /i, u here). Using the form H13|) of the perturbation tensor /i^^ 
we derive the following expression for the curvature terms in the wave equation @ , 

'^Ra^.uphP'' + R^,h\ + R,„h\ = {2a-^a + Aa-^a'' + 6Ka-^)fQ^, . (17) 
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For the covariant d'Alembertian of the perturbation tensor we can write 
in which 



(18) 



Qiiu;oa = (6a ^a^ — 2a ^a)Q^i, , Qi^u-^u — Q^j,u\ii + 2a '^a^^uQ^y 



Combining equations (|17|l . (|18|l . and using the Helmholtz equation (|14|l . the wave 
equation ^ is then reduced to the second-order ordinary differential equation for 
the ampHtude of the perturbations f{rj) depending on the conformal time, 



ad a? 
a ar] 



6K 



(19) 



3.1.1 Waves in the de Sitter universe 

As an example we will now solve the equation l|19() for the de Sitter spacetime 
which is a FRW universe with the cosmological constant A > 0. Due to its unique 
symmetry the de Sitter spacetime admits various specific spacelike sections which 
provide all three cases of the FRW models (fTT|l . (fl^ of constant spatial curvature 
K = 0,+l or -1, see e.g. [19]. 

The spatially flat K = form of de Sitter metric is characterized by the expan- 
sion function 

a{v) - , (20) 
V 

where a = a/S/A and we assume r/ < 0. The minus sign in H20|l guarantees that 
a(t) = aoe*'", where t is the comoving synchronous time. This is the standard 
"inflationary" behaviour with the asymptotics 77 ^ as < — *■ oo. Inserting the 
particular form H2U|I of the function a{r]) into equation (|19|) we arrive at the following 
solution for the time evolution of the waves 



fiv) 



1 



Ci(fc)j(*^,fcry) +C2(fc)Y(z^,fcr;) 



(21) 



where Ci,C2 are arbitrary functions of k, and cylindric Bessel functions of the 
first kind J(i^, z) and of the second kind Y(i^, z) have an imaginary index v = i-^'^. 
This expression is in a full agreement with the result obtained previously using 
synchronous coordinates for the flat de Sitter metric [20]. A typical plot, for fc = 1, 
of the basic modes is given in figure ^ Using the relation between the conformal 
and synchronous times, |7y| ^ e~*/", we can determine the asymptotic behaviour of 
the solution as f — > 00. From the approximate expressions for the Bessel functions, 
J(^, z) ~ z'^ and Y(z^, z) ~ z"'' for z ^ 0, we infer that with i — > cx) the asymptotic 



behaviour is periodic, J 



and Y ■ 



where the specific frequency lu 



|A is the same for all the modes k, depending on the cosmological constant 
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A only. The overall amplitude "envelope" of the waves 1)21(1 grows as e*/^" which 
should be compared to the growth of the background metric functions e^*/". In 
accordance with the cosmic no-hair conjecture the high-frequency waves are thus 
relatively damped. 




Fig. 1. The figure (a) is a plot of the function Re 



is a plot of Re 
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J{i^^,r]) J>, and the figure (b) 



'15 



2 representing evolution of high-frequency gravitational 

waves in fiat de Sitter universe. 



The value K = 1 for metric Hll() . I|12() corresponds to closed spatial sections 
which are three-spheres S^. The expansion function then takes the form 



0(77) 



smyy 



(22) 



and the coordinates cover the whole dc Sitter manifold. Solution to the equation 
p9(l with the expansion function 1(22(1 is 



fiv) 



Ci(fc)pfV3Tfc2-i 

C2(fc)Q(V3Tr ^ 



15 

, I ; cos 77 

2' 2 ' ' 



(23) 



15 



; COST] 



where ^{1/, /i; z) and Q(j^, z) are associated Legendre functions of the first and of 
the second kind, respectively, which satisfy the differential equation 



2 

{l-z^)y" -2zy' +(v{v + \)-^^)y = , 



(24) 



for y{z). The points z — 1^ — l,oo are singularities of this equation (except in some 
special cases) and ordinary branch points of the Legendre functions in the complex 
domain. When we take the branch cuts to be (— 00,— 1) and (l,oo), and if we 
compose the Legendre functions with the cosine function, as in ((23(1 . the singularities 
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occur for 77 = 0, tt. The sample plots of the basic modes of the solution (for fc = 4) is 
presented on figure|21 Using the transformation between the conformal time and the 
synchronous time 77 = 2 arctan(e*/") we may determine the asymptotic behaviour 
of the solution as t 00 which corresponds to rj —^ From the approximations of 
associated Legendre functions as z ^ — 1, namely P(7^, fi; z) ~ Aiy, + z)^/^ + 
B{v, ^)(l + z)-^/2 and Q(i^, ^; z) - C(y, ^^){l + z)^'/^ + D{v, fi){l + z)-f/^, we obtain 
that the asymptotic behaviour of both P and Q is periodic in the synchronous time 

t: it is of the form e^*"* and e*"* as t — > 00, where the unique frequency oj = \J |A 

is again common to all the modes. The amplitude of the waves (|23|l grows as 
^Jcosh{t/a) whereas the de Sitter metric functions grow as cosh^{t/a) so that the 
waves are asymptotically damped with respect to the background. 




Fig. 2. The figure (a) is a plot of the function Re 



(b) 

P(VT9 - \,i^;cosrf) 



and 



the figure (b) is a plot of similar expression in which P is replaced by Q. 



Finally, when K = —\ the spatial sections of the de Sitter spacetime are hyper 
bolic, and 

a-ii) = r 



sinh r/ 



(25) 



Solution to the equation (fl^ with (P5|l is 

1 



sinh I ?7 

+ C2(fc) 



Ci(fc)p(V3-P-i 



; cosh 77 



; cosh 77 



(26) 



Taking the branch cuts again to be (— 00,— 1) and (l,cx)), and composing the 
Legendre functions with hyperbolic cosine, as in (|26|l . we obtain the sample plot 
(for k = 4) shown in figure |3| Considering the relation to a synchronous time, 
77 = In I tanh(<:/2a)|, we observe that the asymptotic behaviour as i ^ 00 corre- 
sponds to 77 ^ 0. If we use the approximation of the Legendre functions P(j^, ^; z) ~ 
E{v,^l){l- zY'"^, q{v,fM;z) ~ F(7/,/i)(l-z)^/2 + G(zy,/i)(l-z)-'^/2 which is valid 
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near z = 1 we again obtain that the asymptotic behaviour of P and Q in the syn- 
chronous time t is asymptoticahy periodic with the common frequency u) = \ |A. 



The amphtude of the waves grows as sinh(i/a) which is slower than the growth 
of the background metric sinh^(t/a). 




Fig. 3. The figure (a) is a plot of the function Re 



^sinh |>7| ^ 2 ' 2 

and the figure (b) is a plot of an analogous expression, only with P replaced by Q. 



:P(V^- |,i^;coshr7) 



3.1.2 Waves in the anti-de Sitter universe 

Analogously, we obtain the evolution of high-frequency gravitational waves in 
the anti-de Sitter spacetime represented by the FRW metric with a negative cos- 
mological constant A. The spacetime can be written in the form (|12|l with 

K = —1 and the expansion parameter 

a{v) = , (27) 

cosh?7 



where /3 — However, these coordinates cover only part of the complete 

manifold. Solution to the equation 1)19(1 with the expansion function H27|l is 



fiv) ^ 



Ci(fc)p(yr~F~i,z^;isinh77) (28) 
+C2{k)q(^V3-k^- i z^;zsinhi 



-^cosh rj 



which is the analogue of (|26|l . The sample plot of the basic modes (for fc = 4) is 
presented on figure^l The maximum of the envelope is obviously achieved at = 0, 
and the waves are damped as |?7| — > oo. 

Due to the specific geometry of the anti-de Sitter spacetime it is not possible 
to put it into the standard FRW form ((11(1 . p2(l with spatial sections K = and 
K = 1. To obtain the evolution of possible high-frequency waves in other standard 
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coordinate forms of the anti-de Sitter universe, one has to start with the general 
wave equation @. Let us illustrate the procedure in the case of the conformally 
flat metric 

ds^ = ^ (-dry^ + dx^ + dy^ + dz^) , (29) 

(these coordinates cover the whole manifold). It is easily seen that using the formal 
complex transformation x — ij], f) — ix, a — i(3 we obtain the conformally flat 
metric (|20() of de Sitter spacetime. This offers the possibility to adopt the results 
obtained for de Sitter spacetime {K = 0), and to find the high-frequency perturba- 
tions for the anti-de Sitter spacetime. Unfortunately, this would imply setting the 
components /ii^ of the perturbation tensor to zero due to gauge condition h^f^ — 
applied in new coordinates. The covariant form of this condition is h^^^v'^ = 0, 
where is the four-velocity of an observer. Therefore, the condition hif^ = 
means that the observer moves faster than the speed of light in the direction of ^ , 
in the coordinates of metric H29|) . 

Hence we will attempt to solve the problem directly using the metric (I29II and 
assuming /i^q = 0. The gauge condition /i^^ = simplifies to /i% — (using the 
fiat space summation convention over the indices of an opposite character). Non- 
trivial components of the gauge condition h^^''^ are xh^j^i — 2hij = 0. Applying 
these conditions on the dynamical wave equations @, for perturbations h^^, we 
obtain 

ft-ll = , /l22 = --^33 , 
^21,0 = , /l3i,o = , ft,21,2 + ft-31,3 = , 

cc" /i2i,fc'= + 4/121 = , a;2 /i3i,fe'= -f 4/131 = , (30) 
2^^(-/i22,oo + /i22,fc'') + 2a;/i22,i - 4a;/i2i,2 + 4/i22 = , 
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a; (-/123.00 + h23,k ) + 2x/l23,l - 2a;(/l3i,2 + /l21,3) + 4/l23 = . 

It is obvious that there are only two dynamical degrees of freedom corresponding 
to /i22 — — /i33 and ft.23- The residual components /i2i and are independent of 
7] and therefore play the role of (supplementary) boundary conditions. The most 
natural choice is to put ^,21 = = h^i. The solution of the set of equations (|30|l can 
thus be interpreted as a purely transversal wave which propagates in the direction 
of ^ and which has two polarizations. By this method we do not obtain waves 
propagating in a general direction — this is a consequence of "anisotropy" of the 
anti-de Sitter spacetime in the coordinates H29|) . 

With the above choice, the unified form of wave equation for both degrees of 
freedom /122 = —^33 and /123 is 

drf dx^ dy^ dz^ x dx ' 

where / stands for /122 or /i23. Considering the separation of variables / = g{x) 
X exp[j(— fcof? + + ^3^)] we obtain 

ff" + ^5'+(^ + fc?)3 = 0, (31) 

where kl — fep — — Equation H31() is formally equivalent to H19() for the de Sitter 
spacetime with K = Q (when replacing x with 77, and k\ with fc^, so that a — (3/x). 
The general solution is thus a monochromatic high-frequency gravitational wave of 
the form 

exp [ii-koT] + k2y + k^z)] . 



Ci(A:i) fcix) +C2(fci) y(i^, fcia;) 



(32) 

Let us finally mention that Isaacson's high-frequency waves in the anti-de Sitter 
spacetime can also have the form of the Defrise family of exact solutions [18,21]. 
These metrics were investigated in [22] using the coordinates 

ds^ = /32(d6'2 + sinh^ ed<j)^) + 8/32(cosh6l + sinh 6I cos 0)2dudw 

-Hl6/32(cosh6i + sinh6lcos(/))^d(u)du2 . (33) 

The wavefronts u — const, are two-dimensional hyperbolic surfaces with con- 
stant negative curvature —f3 parameterized by 9 and 0. The anti-de Sitter back- 
ground is given by d(u) = 0, and since the component 7„„ of the metric is of the 
form R{9, (j))d{u)^ it may also represent high-frequency perturbations for small but 
rapidly varying function d{u) = 0{e). The gauge conditions ([7|l, © are fulfilled 
identically while the wave equation @ is satisfied to the order 0(e). 

3.2 Waves in the anisotropic Kasner universe 

We also investigate Isaacson's high-frequency gravitational waves in cosmological 
model with lower symmetry, namely in the Kasner universe. This is a special case 
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of the Bianchi type I class of homogeneous but anisotropic spacetimes. Its metric 
in synchronous coordinates has the standard form [18] 

ds^ = -dt^ + t^'P'dx^ + t^PMy^ + t^P^dz^ , (34) 

where pi,p2,P3 are constants, which represents a solution of vacuum Einstein's 
equations if 

Pl+P2+P3^l, Pi^ +P2^ +P3^ = 1 ■ (35) 

However, we need not assume these relations. It is possible to consider the matter 
content of the universe described by the energy-momentum tensor that does not 
contain a derivative of the metric. This would fulfill the conditions of a generaliza- 
tion of the Isaacson approximation method to non- vacuum spacetimes, as described 
in [13]. Let us mention that it has recently been shown [23,24] that it is impossible 
to retain anisotropy when the Kasner universe is filled with a viscous fluid, such 
that dominant energy condition holds and entropy is non-decreasing. Nevertheless, 
anisotropy is permitted when it is filled with an ideal fluid satisfying the Zel'dovic 
equation of state. 

As in the previous calculations we will use h^o = as an additional condi- 
tion. The traceless gauge condition (jSJ has the form (using the spatial summation 
convention) 

t~^P^h\ = . (36) 
The gauge condition |7l gives the equations 

p^t-^P'h\ ^ , r'^P^h'j^, = . (37) 

By applying these conditions we simplify the dynamical equations © to 

Pit-'^P^h'j^, = and 

- %,oo + {-"^Pk + 2k + 2pj^th^j^o + t^'^'^-P^'^h^j^k'" - ■ip.pj h,j = (38) 

k 

(no summation over i,j here). Using the covariant d'Alembert operator, the differ- 
ential equation (|38|l can be rewritten as □ hij + 2(pi+pj) t~^hijX) — '^PiPj t~^hij = 0. 
Let us denote an arbitrary component hij of the perturbation tensor simply as / 
(even though the wave equation is different for different indices and let us 
define 

A^2{p,+pj)-^Pk , p,p, . (39) 

k 

We will look for the solutions of (|38|l in the following special form, 

f{t, X, y, z) = X(t, x) -f y(t, y) + Z(t, z) . (40) 

Inserting H4()|l into the wave equation H38|l . its left-hand side splits into three parts, 
each depending only on one spatial coordinate. The simplest possible solution is to 
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equate each of these parts to zero, satisfying thus the equation. For example, in the 
spatial coordinate x we obtain 

d'^X AdX ^2 



AB 



■X ^0 



(41) 



9*2 t dt dx^ 
Applying now the one-dimensional Fourier transform in the coordinate x on the 
equation (|41|l we arrive at the ordinary differential equation 



dt^ t dt +^ k,X+ ^^X-0 



(42) 



where X = J^[X]. Making an ansatz X = ti''^+'^^ F{t) and using the coordinate 
transformation s — ji^fci t^^^^ we obtain the standard form of the Bessel equation 

165- (A+ 1)2 



s'G + sG 



G = , 



(43) 



where G(s) = F{t), and the dot denotes differentiation with respect to s. Solving 
H43|l and transforming this back to X and t, the solution of evolution equation H42(l 
takes the following form 



i+(fcl)j( 

+cr(fci)Y( 



2(1-Pi) 



2(1-Pi) 



l-pi 



(44) 



For vacuum Kasner universe (for which the relations H35(l hold) we obtain (j4+1)2 — 
165 = 4(pi — Pj)2, and thus the index of the Bessel functions is a real number. 
Generally, assuming that all pi are positive, it turns out that for J^k Pk < ^ the 
index is always real, but for > 1 it may be imaginary. 

Proceeding in the same way for the functions Y and Z, the form of the solution 
is reproduced except for the replacement of fci,pi, C^, Cj" with k2,p2, Cj^, , 
and k^,p3,C^ , respectively. The complete solution to H38|) may thus be com- 
posed in the following way 

fit, X, y, z) = [Xit, ki) S{k2)S{k3) + Y{t, k^) S{ki)S{k3) + Z{t, k^) 5{k^)6{k2)\ , 

(45) 

where J-3 ^ denotes the inverse Fourier transform in three dimensions, and 5 de- 
notes the Dirac delta function. The high-frequency gravitational waves are then 
determined by the six functions C^{kj), j = 1, 2, 3. The monochromatic wave with 
the wave- vector (fc", ^3) obtained by setting Cf{kj) = cf5{kj — fc"), and has 
the form 



/(i,xi,X2,a;3) = i^(^+^) ^ c+j( 



2(1-P3) 



1 



(46) 



V (A+l)^-16a k]t^^-P^') 



exp {ik'^Xjj 
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where xi, X2,X3 stands for x, y, z. 

Since the Bessel function Y is generaUy divergent when approaching zero, there 
obviously arises the question of consistency of the above perturbative approach 
at vicinity of i = when pi < 1. Asymptotic behaviour of this Bcssel function is 
Y{i',z) ^ for z ^ 1. Analysis for general Kasner solutions is difficult due to 
the complicated dependence on three parameters . However, two particular cases 
can easily be investigated: 

— For purely imaginary index v = ia\t behaves as a complex unit and therefore 
it is non-divergent: 

Y(y,z) - e-'"'"^ . 

— For vacuum Kasner universe = 1) one calculates that 

X ^ i3[^+l-\/('4+l)"-16B] ^ l{{v^^V3)-\v^-V^\ = ^Zp ^ 

where p = min{pi,pj} , 

which corresponds to the behaviour of the metric coefficients ~ t^'^i in H34|) . 
The Isaacson approximation is thus consistent since the amplitudes of the 
perturbations (i.e., c^tCJ) are small. 

4 Conclusions 

We investigated time evolution of Isaacson's high-frequency gravitational waves 
in some standard cosmological models. We concentrated mainly on homogeneous 
and isotropic FRW universes by decomposing the perturbations into tensor har- 
monics. The Isaacson approach takes the influence of high-frequency waves on the 
backgrounds into account which is the main difference to other studies of gravita- 
tional waves in the above cosmological spacetimes. 

Illustrative explicit results were presented for the maximally symmetric de Sitter 
and anti-de Sitter spacetimes in various forms having different spatial curvatures 
K — 0^ 1,-1. The complete solution representing evolution of gravitational waves 
which propagate in an arbitrary direction was obtained — in terms of the Bessel 
or Legendre functions — for spacetime metrics with isotropic time-slices. Interest- 
ingly, since the constant parameters a and (3 occur in the conformal factors 0(77) 
given by (^0)) , , (^5)) , l(?7|) only multiplicatively, the corresponding equation ((T^ 
expressed in the conformal time 77 is independent of the cosmological constant in 
these cases. This explains why the resulting spectra H21II . (j22J, l|26|l . (|28|l do not 
explicitly depend on the value of A. However, if one introduces the corresponding 
synchronous time, the dependence on the cosmological constant occurs. For exam- 
ple, in the de Sitter models the Isaacson wave perturbations for all modes k and 
for any spatial curvature K of the background become periodic asymptotically as 

t ^ 00 with the common characteristic frequency uj — ^ |A. Moreover, the waves 
are exponentially damped — relative to the expanding de Sitter background — in 
full accordance with the cosmic no-hair conjecture. 
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High-frequency gravitational waves in the anisotropic Kasner universe were also 
studied assuming perturbations with a decoupled spatial dependence. In this case 
the possible spectrum can be expressed in terms of the Bessel functions. We an- 
alyzed the evolution of these perturbations and demonstrated (assuming pi < 1) 
that near the initial singularity t = the Isaacson approximation method is consis- 
tent. For t — > oo the vacuum perturbations, relative to the background, disappear 
so that the spacetimes evolve to homogeneous ones. 
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